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Abstract 

We include alignment interactions in a well-studied first-order attractive-repulsive 
macroscopic model for aggregation. The distinctive feature of the extended model is 
that the equation that specifies the velocity in terms of the population density, becomes 
implicit, and can have non-unique solutions. We investigate the well-posedness of the 
model and show rigorously how it can be obtained as a macroscopic limit of a second- 
order kinetic equation. We work within the space of probability measures with compact 
support and use mass transportation ideas and the characteristic method as essential tools 
in the analysis. A discretization procedure that parallels the analysis is formulated and 
implemented numerically in one and two dimensions. 

Keywords: aggregation models; nonlocal interactions; kinetic equations; macroscopic 
limit; mass transport; particle methods 


1 Introduction 

The literature on self-organizing behaviour or swarming has grown dramatically over the last 
years. A variety of mathematical models has been proposed, which have origin in biologi¬ 
cal applications (e.g., self-collective behaviour seen in species such as fish, birds or insects 
[CDF'*'03]), as well as in social sciences and engineering (e.g., opinion formation [MT14], so¬ 
cial networks [JaclO], robotics and space missions [JE07]). The main aspect is the modelling 
of the social interactions between the members of a group; due to such inter-individual inter¬ 
actions, self-organization may occur in a physical space (insect swarms, fish schools, robots) 
or, more abstractly, in an opinion space. 

One approach in modelling aggregation is to consider individuals/organisms as point par¬ 
ticles and design either an ordinary differential equation (ODE) or a discrete-time equation 
to model their evolution. Another is to formulate a partial differential equation (PDE) that 
governs the time evolution of the population density field. These two approaches result 
from the various descriptions that one can take in modelling aggregation behaviour: particle- 
based/microscopic or continuum/macroscopic. We refer to [CFTVIO] for a recent review of 
aggregation models, where in particular, it is shown how microscopic models can be related 
to macroscopic ones via kinetic theory. 
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Three types of social interactions have been commonly considered in the literature on 
mathematical aggregations: attraction, repulsion, and alignment. Consequently, aggrega¬ 
tion models can be distinguished in terms of which of these interactions are being accounted 
for. Some models consider only a subset of these interactions (just attraction and repulsion 
[MEK99, BCL09, DCBC06, BTll] or just alignment [CS07, CFRTIO]), others account for all 
three of them. Models of the latter type are typically referred to as “three-zone” models, 
as each particular interaction type acts at different ranges (repulsion acts at short distances, 
while alignment and attraction are present at intermediate and long ranges, respectively). 
This class of models has had many successful applications in biological and sociological mod¬ 
elling [CKJ+02,Rey87,MKJNR12]. 

Aggregation models (discrete or continuous) may also differ in how the velocity field is 
determined. There are second-order/dynamic models, typically in the form of the Newton’s 
second law, where a differential equation for the evolution of the velocity is being provided 
[DCBC06, CS07], and first-order/kinematic models where the velocity is prescribed in terms 
of the spatial configuration [MEK99, TB04, BCL09]. The aim of the present paper is to 
extend, by including alignment interactions, a first-order continuum model for aggregation 
that attracted a high amount of interest in recent literature [MEK99, BV06, BL07, KSUBll, 
LTB09, FHKll]. Below we introduce the extended model and its derivation, then point out 
the fundamental issues that arise with such an extension, and how we address these issues in 
the present paper. 

Consider the following continuum model for the evolution of the macroscopic density 
function p{t,x) in R'’*: 

dtp + • {pu) = 0, = po{x), (1.1a) 

d>(t, x)u{t, x)= (/>(|x - y\)p{t, y)u{t, y) dy - * p{t, x), (1.1b) 

R'i 

where (j) is an influence function that controls the alignment interactions, ^ = (p*p, and K is an 
attractive-repulsive interaction potential. The asterisk denotes spatial convolution. Hence, 
the model consists in an active transport equation for the density p, with velocity field u 
defined by (1.1b). The coefficient <I>(t, x) in the left-hand-side of (1.1b) has the interpretation 
of the total influence received at location x and time t from the rest of the group. The 
right-hand-side of (1.1b) has two terms: the first models alignment, and the second models 
attraction and repulsion, as detailed below. 

Alignment is modelled through an averaging mechanism that allows individuals to adjust 
their velocities relative to the velocities of the others. Specifically, the velocity at location 
x and time t is assumed to depend non-locally on the velocities u{t, y) at locations y within 
the alignment interaction range set by the support of the influence function (f)- The first-term 
in the right-hand-side of (1.1b) captures this averaging process, with weights/interaction 
strengths given by (j){\x — y\), assumed to depend only on the relative distance between 
locations x and y. 

Attraction and repulsion are modelled by the convolution of the gradient of the interaction 
potential with the population density. In brief, individuals are assumed to repel each other at 
short ranges, to create a comfort zone around them, but attract each other once they distance 
themselves too far apart. Equation (1.1a) with the velocity field given solely by this term, i.e., 
u{t,x) = —VxK * p{t,x), constitutes the aggregation model that has been referred to above 
and which the present paper generalizes. A variety of issues has been investigated during 
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the last decade for this explicit attractive-repulsive aggregation model, including the well- 
posedness of solutions [BV06, BL07, BLRll, CDF+ll], the long-time behaviour of solutions 
[MEK99,BDF08,KSUB11,LTB09,BT11,FHK11], and the derivation of the continuum model 
as a mean-field limit [CCH14]. 

We point out an essential feature of equation (1.1b), which is that it is an implicit equation 
in u, and can have non-unique solutions (e.g., due to translational invariance, one can add an 
arbitrary function of t to any solution of (1.1b), and obtain a different solution). This is a key 
challenge brought up by the inclusion of alignment interactions in the explicit aggregation 
model from [MEK99]. Addressing this challenge is one of the major goals of the present paper. 
The non-uniqueness of solutions to models of type (1.1) has been noted in [MKJNR12], but 
no resolution was offered. To our best knowledge, the present paper is the first systematic 
study of a first-order continuum model for aggregation that includes both attractive/repulsive 
and alignment interactions. 

The origin of the macroscopic model (1.1) can be traced back to the following second-order 
discrete model derived from Newton’s second law. Suppose there are N particles in M'’*, whose 
positions and velocities denoted by Xi and Vi, respectively {i = 1,, N), evolve according to 
the following system of ODE’s: 

dxi 

for i = 1,... ,N. Here, it has been assumed that that all particles have the same mass rrii = e. 
The functions and K have similar meaning as in (1.1). 

Without the attractive-repulsive term modelled by the second term in the right-hand- 
side of (1.2b), model (1.2) represents the celebrated model of Cucker and Smale [CS07]. 
It is well-known that for certain influence functions (j), the Cucker-Smale model successfully 
captures the unconditional flocking phenomenon, where individuals align their velocities into a 
certain asymptotic direction [HT08,HL09,CFRT10]. Comprehensive literature also exists on 
second-order attractive-repulsive models without alignment [DCBC06, CDM'’'07, CETVIO]. 
Second-order models with both alignment and attraction/repulsion have also been studied 
[CD10,LLEK08,CCR11, AIRll], though in not as much depth and detail as the models with 
the two sets of forces considered separately. 

The passage from the second-order discrete model (1.2) to the first-order macroscopic 
model (1.1) is done in two steps. First, for each fixed e > 0, one can take the limit —)• 00 in 
(1.2) and reach, by BBGKY hierarchies or mean field limits (see e.g., [HT08] and the review 
in [CFTVIO]), the following kinetic equation for the density f{t,x,v) at position x G and 
velocity u G 

dtfe + v-V^fi = ^Vy-{F[fi]fl), fe\t-^Q = fo{x,v), (1.3a) 

^[fe] = J 4>{\x-y\){v-v*)fl{t,y,v*)dydv* + J VxK{x - y)fi{t,y,v*)dydv*. (1.3b) 

R2d ]R2d 

Rigorous derivations of mean field limits starting from particle systems is a classical subject, 
comprising an extensive body of works. Some of the most recent works include the mean field 


(1.2a) 

(1.2b) 
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limit of the Cucker-Smale model [HL09, CFRTIO], as well as extensions of these results to 
include general aggregation models of the form (1.2) [CCRll]. 

The second step in passing from (1.2) to (1.1) is to send e ^ 0 in the kinetic equation (1.3) 
and derive (1.1) as a hydrodynamic limit. The rigorous treatment of this limit constitutes in 
fact a major component of the present work. Passage from kinetic to macroscopic equations 
is also a vast topic, extensively studied for instance in the context of hydrodynamic limits 
of the nonlinear Botlzmann equations. It is beyond the scope of this introduction to give a 
detailed account of this well-established research area, we simply refer here to a recent review 
paper [SR14] and the references therein. 

As indicated above, a major issue that arises when one considers the first-order model 
(1.1) is the non-uniqueness of solutions to (1.1b). In the present paper we resolve this non¬ 
uniqueness issue for the case when the interaction potential is symmetrie about the origin, 
i.e., it satisfies 


K{x) = K{—x), for all x E 


(1.4) 


For symmetric potentials, the ODE system (1.2) and the kinetic equation (1.3) conserve the 
linear momentum: f^ 2 dvfe{t,x,v)dxdv, respectively. In terms of macroscopic 

variables, J^dPeUedx remains constant through the evolution of (1.3), where and Ue are 
dehned by: 

Pe{t,x)= j fe{t,x,v)dv, p^it, x)u^{t, x) = J V f^it, X, v) dv . 

The key idea for dealing with the non-uniqueness of solutions to (1.1b) is to account for the 
fact that (1.1) is the limit of a second-order model for which linear momentum conservation 
holds and hence, as a limiting equation, (1.1) should inherit this momentum conservation as 
well. 

Given the considerations above, we append model (1.1) with the constraint/requirement 
that the linear momentum remains constant through the time evolution. With no loss of 
generality, we assume that the linear momentum is zero at the initial time and hence, it has 
to remain zero for all times. For ease of referencing, we list below the macroscopic model 
(1.1) together with momentum conservation: 


dtp + • {pu) = 0, = po{x), 

4>(t, x)u{t, x) = J (pilx - y\)p{t, y)u{t, y) dy - * p(t, x) 


/ 


pudx = 0. 


(1.5a) 

(1.5b) 

(1.5c) 


System (1.5) is the main object of study of this paper. 

The layout of the present work is as follows. In Section 2 we establish a well-posedness 
theory for system (1.5). Solutions are sought in the space of probability measures 

with compact support, using the mass transportation framework developed in [CCRll]. The 
main difficulty is to obtain a Lipschitz bound on u, given that for each p, (1.5b) provides 
multiple solutions. The momentum conservation condition (1.5c) plays an essential role, as 
it enforces uniqueness. The well-posedness result is stated in Theorem 2.4. 

Another major issue of interest investigated in detail in this paper is the zero-inertia 
limit of solutions of the kinetic model (1.3). We use again the framework from [CCRll] 
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and consider measure-valued solutions of (1.3) in the mass transportation sense. For this 
reason, the method of characteristics plays a major role throughout our analysis. In Section 3 
we establish uniform bounds in e of solutions to (1.3). We note in particular the key bound on 
the second-order fluctuation term established by Theorem 3.7. The uniform bounds enable 
us to derive the main convergence results in Section 4. Specifically, we show that converges 
weak-* as measures to f{t,x,v) = p{t,x)5{v — u{t,x)) as e —)• 0, where p (and u) represent 
the unique solution to the macroscopic model (1.5). The precise statements are presented in 
Theorems 4.1 and 4.5. We also point out the convergence of characteristics paths established 
by Theorem 4.4, where the subtlety resides in the singular limit, with characteristics of a 
second-order system collapsing into first-order characteristics. 

Finally, in Section 5 we design and implement a numerical procedure for solving (1.1). The 
numerical method is rooted in the analytical considerations made for the continuum model, 
as a discrete analogue of the momentum conservation is enforced in order to obtain unique 
solutions for the discretization of (1.5b). To test our method, we implement and illustrate it 
in one and two dimensions, for several choices of interaction kernels K for which we know the 
exact steady states. 

We note that for models that include only attraction and repulsion, the analogous zero 
inertia limits e —)• 0 of second-order models^ such as (1.2) and (1.3), have been investigated 
recently in [FS14] (an earlier study, restricted to PDF analysis, was done in [JabOO]). As in 
[FS14], the focus in the present work is to find uniform (in e) bounds on the size of the support 
of the solutions /e, as well as to show the vanishing (with e) of a second-order fluctuation 
term. By adding the alignment term such estimates become significantly more involved, as 
they require the use of the momentum conservation in (1.3), as well as the coupling with 
control on high order moments of /g. 

Finally, we point out that we work in this paper with smooth interaction potentials and 
influence functions. Specifically, we assume ViF G VF^’°°(M'^) and (j) € VF^’°°(M^), where (j) is 
also assumed to be positive, non-increasing, and to decay sufficiently slow at infinity. All these 
assumptions are used in an essential way in our analysis, in particular to control (uniformly 
in e) the support of the solutions to (1.3) and the second-order fluctuation term. Removing 
(some of) these requirements would require a significantly different approach and analysis 
than the present one. We comment however that, from the point of view of applications, 
the smoothness requirement is mostly irrelevant, as in numerical simulations one does not 
normally distinguish between a singular potential and its smooth regularization. 

2 Well-posedness of the macroscopic equation 

In this section, we establish a well-posedness theory for the macroscopic equation (1.5). Solu¬ 
tions p are sought in the space C([0, T]-,Vc{^'^)), where Rc(I^'^) denotes the space of probability 
measures with compact support in R'^. We use the measure-transportation framework from 
[CCRll] and consider the characteristic equations associated to (1.5a), given by 

-^x = u[p]{t,x{t)), x{0)=xo. (2.1) 

In order to specify the solution space for (1.1), we use the concept of flow maps: suppose p is 
a fixed function and the vector field u[p\{t,x) is locally Lipschitz in x. Then standard ODE 

^We refer to the kinetic model (1.3) as being “second-order” since it is based on Newton’s second law (1.2); 
in strict terms, equation (1.3) is a first-order PDE. 
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theory provides the local existence of the flow map associated to (2.1): 



where x{t) is the unique solution of (2.1) that starts at xq. 

Definition 2.1. Let T > 0 be fixed. A measure-valued function p G C([0, T]; Pc(lK'^)) is said 



( 2 . 2 ) 


where at each t G [0,T], the velocity field u[p]{t, ■) G W^’°°{suppp) satisfies (1.5b) and (1.5c). 


The push-forward in Definition 2.1 is taken in the mass transportation sense; equation 
(2.2) is equivalent to: 



for all If G C'fe(M'^). 

The results in [CCRll] establish the existence and uniqueness of measure solutions for 
various aggregation models via fixed point arguments. These results apply in particular to 
attractive-repulsive models in the form (1.5a), but there the velocity field u is given explicitly 
in terms of p. In [CCRll], the assumptions made on the attractive-repulsive potential K 
guarantee the local Lipschitz continuity of the velocity field u and hence, the global existence 
and uniqueness of characteristics paths. 

In our model, u is given implicitly by (1.5b) and (1.5c). While (1.5b) can have multiple 
solutions (as discussed above), we show that the momentum conservation property (1.5c) 
enforces uniqueness. Once a unique solution u[p\ of (1.5b)-(1.5c) has been identified, this 
solution is shown to be Lipschitz continuous in the spatial variable so that the characteristic 
equation is well-posed. 

We now establish the necessary a priori bounds for (p, u) for deriving the well-posedness 
of solutions to (1.5). First we show the Lipschitz bound of u for a fixed density function p. 
Similar estimates have also been used in [MTll,TT14] for flocking models. 

Proposition 2.1. Assume K is symmetric (i.e., it satisfies {l.A)), S/xK G 1F^’°°(M'^) and 
4> G VF^’°°(M"'') is positive and non-increasing. Consider a given probability measure p{t, •) G 
Rc(l^^) (it some fixed time t. Then equations (1.5b)-(1.5c) uniquely define a bounded, Lipschitz 
continuous velocity field u on the support of p. 

Proof. Since equations (1.5b)-(1.5c) hold at a specific fixed time, we simplify the notation, 
and drop the time dependence in the proof. Denote 


Sp = max \x\. 


xGsuppp 


(2.3) 


By our assumption on p, Sp is finite. The following lower bound on the influence function 
can be inferred immediately: 


fi{\x — y\)>ri for any x,y G suppp. 
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where 


(2.4) 


rj := (^(25p) > 0. 

As a consequence, <h(x) is bounded from below. Indeed, for any x E supp/o, 

4>(x) = J 4>{\x - y\)p{y)dy >r]. (2.5) 

suppp 


Moreover, we have 

^(a;) < \\4>\\ L°°(o,oo) > foi' any X E 
Define the following alignment operator A: 

A[u]{x) = <l>{x)u{x) - J 4i{\x -y\)p{y)u{y)dy. (2.6) 


It is easy to check that M is a linear bounded operator which maps L°°(suppp) to (suppp). 
Equation (1.5b) can now be expressed as 

M[m](x) = -VxK * p{x). 

As commented above, the operator A is not invertible since A maps all constant (in x) 
functions to zero. 

To overcome the degeneracy of M, we make use of the momentum conservation condition 
(1.5c) and define a new operator A4: 


M[u]{x) : = A[u]{x) + p / p{y)u{y)dy 


^{x)u{x) - J (())(|x -y\)- v)piy)u{y) dy. 


(2.7) 


We note that the operators A and Ai depend on p; however, we choose not to indicate this 
dependence explicitly in the notations, as this fact is irrelevant for the main well-posedness 
result that follows. 

A solution u of (1.5b) and (1.5c) also satisfies 

Ai[u]{x) = —VxK * p(x). (2.8) 


We claim that Ai is invertible. Indeed, given any x E suppp, 


J i4>i\x - y\) - rj)p{y)u{y)dy < \\u\\l^^suppp) j {4>i\x - y\) - v)piy)dy 

Rd suppp 

II^IIl°°(suppp) (^(®) V )' 

By the definition of Ai in (2.7), one has 

Ai[u\{x) > ^{x)u{x) — ||u|| 1^00 (suppp) - rf) for all x E suppp. 
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For any <5 > 0, there exists G suppp such that u{xs) > ||u|| 2 ,<x>(suppp) — 5 and hence, 

[it] (a;,5) > ^(^< 5 ) ||'1^||l°°(suppp) II'*^IIl°°(suppp)(^(^5) ^ ll*i^llL°°(0,co) 

= V II'*^IIl°°(suppp) ~ ^ II‘?^IIl°°(0,oo) • 

By letting (5 —>• 0 we obtain 


11-'^ [^] ll-I/°°(suppp) — ^11 (suppp) • (^•9) 

Hence for a given p, A4 is invertible on L°°(suppp) and u is uniquely solvable from (2.8) such 
that 

u{x) = {u[p\){x) = -M~^ iVxK * p{x)). (2.10) 

Next we show that u G VF^’°°(supp/9) and find an explicit bound on its VF^’°°-norm. We 
start with the L°° estimate. From (2.9) and (2.10), we have 

II'“IIl°°(suppp) — 11-^ II ll^a:-^ * P|Il°°(suppp) — ~ || Vr . (2-11) 

Next, we estimate the bound of VxU. The distributional derivative of u satisfies 

v,u = piyMy)dy + ^x ■ 


Applying the lower bound estimate (2.5) for and the L°° estimate (2.11) for u in each term 
gives 




(t>{\x - ?/|) 

<l>(x) 


p{y)u{y)dy 


< 


2 nv, 


— 9 liv- ||L°° ll VcA'Iilo 

T]^ 


(suppp) 






^ ~ll'/’^lk°°l 

L°° (suppp) y 


1 

y 


'^xK\\l°° h—II^i-^IIl°°- 


This shows 


3 1 

|VxU||ioo(guppp) < '^\\ 4 > ||l°° II Va;iF||L°° + “l|yr-^l|L°°- 


Combining (2.11) with (2.12), we have 


u 


<(- + 4ii0' 


livi-°°(suppp) — ' ^2 IIl°°(0,oo) I II IlL' 

with p defined in (2.4). 


VxK\\^_oa 4 -Vj-iF , 


( 2 . 12 ) 

(2.13) 

□ 


Remark 2.2. The values of the velocity field u outside the support of p do not explicitly 
enter in the dynamics of model (1.5). However, we can extend the definition of u outside 
of suppp by using (1.5b). This extension is applied in the proof of Theorem 2-4 where a 
W^'°^-hound on u is used in a domain that is larger than suppp. In particular, we note that 
the right-hand-side of (1.5b) only depends on values of p and u inside the support of p. For 
any R > 0 such that suppp C Br and x G Br \ suppp, <l>(x) is bounded from below by (p(2R) 
(see (2.5) j and the assumptions made on K and cj) along with the bound (2.13) immediately 
imply a bound on that depends on R, || Va;iF||^yi,oo and || 0 ||p[/i,oo This bound has 

a similar expression as in (2.13), but with p replaced by 4>{2R). 














As discussed below, Proposition 2.1 combined with the fixed-point argument developed 
in [CCRll] is enough to conclude the local existence of solutions to (1.5) in the measure- 
transportation sense of Definition 2.1. In order to show the global existence, we need to prove 
a priori that the support of p does not become infinite in a finite time. This requires an 
additional assumption on the influence function (p. The precise statement is as follows: 

Proposition 2.3. Assume K and p satisfy the hypotheses of Proposition 2.1, together with 
the following slow deeay condition on the influence function cp: 

OO 

J (p{r)dr = oo. (2-14) 

0 

Let {p,u) be a solution to (1.5) with eompaetly supported initial density pQ. Then for any 
t > 0, we have 

(a) p{t,-) remains eompaetly supported and the bound of its support satisfies 

Sp{t) < (4'(25p(0)) + 2t ||V,K||^^), (2.15) 

where Sp{t) is defined in (2.3) and 'l'(r) = 0(s)ds. 

(b) u G W^’°°{suppp) and (2.13) holds with 

rj{t) = cP{^-\d>{2Sp{0)) + 2t IIV^ATII^^)). 

Proof. First, from (2.11) we know that the L°°-bound of u depends on Sp. Specifically, 

ll'“(^>-)llL°o(s„ppp(t)) < ^^2Sp{t)) 

for all t in a local interval of existence [0,T]. On the other hand, for any characteristic path 
x{t) originating at xq E supppo (see (2.1)), we have 

— |x(t)| < ll'lt(^) ■) IIl°°( suppp(t)) • 

Hence Sp is bounded in terms of u as 

t 

Spit) < Sp(0) + j ||^^(r,-)llL-(suppp(r)) dr. 

0 

We can now close the loop between Sp and u. To this end, denote 

t 

D{t) = 5p(0) + J ||u(r, •)llL-(suppp(r)) dr. 

0 

Then Spit) < Dfi). By (2.16) and the monotonicity of (p, 

D it) = ||M(t, •)llL°°(suppp(t)) < (p{2Dit)) ’ 
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which implies 


By integrating the equation above, we find 

^(2L>(t)) < ^'(25p(0)) + 2t ||V,,i^||^oc . (2.17) 

Assumption (2.14) implies that 'l'(r) —)• oo as r —>■ oo. Moreover, 'h(O) = 0 and 4’^ > 0. 
Consequently, 4' is a bijection from M'*' to M^. We can thus infer (2.15) from (2.17) and the 
inequality Sp{t) < D{t). The Lipschitz bound on u in part (b) follows immediately. □ 

Combining the a priori bounds in Proposition 2.3 with the fixed-point argument in [CCRll], 
we can now show the global wellposedness of solutions to (1.5). The theorem states 

Theorem 2.4 (Well-posedness of the macroscopic model). Assume K is symmetric, VxK G 
1T^’°°(M'^) and (/> G 1T^’°°(]R''') is positive, non-increasing, and satisfies the slow decay con¬ 
dition (2.14). Suppose the initial measure po G Vc(R'^)- Then for any T > 0, there exists a 
unique solution {p, u) to (1.5) such that p G (^([0, T], Pc(M'^)) and u G <^([0, T]; {suppp)). 

Proof. The proof follows directly from the theory developed in [CCRll]; we only summarize 
the setup and the main steps. We refer for details to the proof of Theorem 3.10 in [CCRll]. 

Local well-posedness. Fix R > 0 such that the ball Br contains the initial support of po, 
and consider the set S consisting of all functions in C([0,T],Pc(lK‘^)) such that the support 
of pt(-) = p{t,-) is contained in the ball Br for all t G [0,T] (T is a positive number to be 
chosen later). The set S endowed with the distance Wi defined by 

y^i{pt,Pt)= sup Wi{pt,pt), 

te[o,T] 

is a complete metric space. Here, Wi denotes the 1-Wasserstein distance. 

For fixed p & S, one can define a map Q on this metric space by 

Consequently, solutions of (1.5) in the sense of Definition 2.1 are fixed points of the map Q. 
The essential step is to show that this map is contractive, and hence, conclude that it has a 
unique fixed point in S. 

In our context, the key observation is that for every p ^ S, one can use the argument 
in the proof of Proposition 2.1 (with p = 4>{2R) independent of t G [0,T] and p G S) and 
find a VF^’°°-bound on the velocity held u[p], bound which depends only on K, and R (see 
estimate (2.13) and Remark 2.2). 

The Lipschitz continuity of u[p] provides well-posedness of characteristic trajectories for 
each p G S. Also, for T short enough, the support of Q[p]{t) is contained in Br. By simple 
arguments (Lemmas 3.11 and 3.7 in [CCRll]), the following estimate can be shown to hold: 

Wi{g[pi]{t),g[p 2 ]it)) <c{t) sup \\u[pi]-u[p 2 ]\\l--(Br)- 

te[o,T] 
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To infer the contractivity of Q we also need an estimate of the form 


sup \\u[pi]-u[p 2 ]\\l^,s )<C{T) sup Wi{pi,p 2 ), (2.18) 

t£[0,T] ' t£[0,T] 

for some constant C that depends on T, as well as on K, cf) and R. Requiring that T is small 
enough such that C <1 concludes the existence of a fixed point of Q, and consequently, of a 
local solution to (1.5) in the sense of Definition 2.1. 

For proving (2.18), let pi, P2 £ <S, and rj = (p{2R). Denote ui = u[pi] and U2 = u[p2\- 
Then, by (1.5b) and (1.5c), 

{(p* pi)ui{t,x) = j- y\) - r])pi{t,y)ui{t,y)dy - * pi{t,x), 

Br 

{ 4 >* P2)u2{t,x) = J { 4 >{\x - y\) - r])p2it,y)u2it,y)dy - * p2{t,x). 

Br 

By Remark 2.2, we can extend ui and U 2 to the whole ball Bn such that ui,U 2 G 
Taking the difference of the above two equations, we have 

{(p * pi){ui -U2) = J {(pi\x - y\) - v)pi{y){ui{y) - U2{y)) dy + 7 ^l, ( 2 . 19 ) 

Br 


where 




-VxK * {pi - P2) 


{(p * {pi - P2))U2 + 



Br 


y\) - v){p\{y) - P2{y))u2{y) dy. 


By Lemma 3.15 in [CCRll], we have 

— (ll^:E-^IL°o(Kd) + ll'?^llrvi.°°(0,oo) Wi{pi, P2) , 

where a bound on \\u2\\\Yi,^(^nR) ^e inferred as explained in Remark 2.2. 

From (2.19), following an argument similar to that leading to (2.11) from (2.7) and (2.8), 
we estimate 


||ui(t, •) -M2(L-)llL°°(Sfl) ^ CWi{pi,p 2 ), for all t € [0,r], 

where C depends on T, R, 00 ) Hence (2.18) holds and there 

exists T > 0 small enough such that (1.1) has a unique solution {p,u) over the time interval 

[o,r]. 

Global well-posedness. Note that the argument above does not use the slow decay condition 
(2.14) on (p. Local existence can be extended for as long as the support of p remains bounded. 
Provided that (p satisfies the additional assumption (2.14), Proposition 2.3 guarantees that 
the support of p does not blow up in a finite time. In this case the local theory can be 
extended to [0, T] for any T > 0, which gives the global well-posedness. □ 
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Remark 2.5. It is interesting to note that the slow decay condition (2.14) appears in previ¬ 
ous literature on the Cucker-Smale model, in connection to asymptotic flocking. Indeed, for 
influence functions in power-law form, it has been shown in [HL09] that the discrete Cucker- 
Smale model has the unconditional flocking property when (j) satisfies (2.14). In addition, 
the same property (2.14) is needed for unconditional flocking in the corresponding kinetic and 
macroscopic Cucker-Smale models [CFRTIO, Tanl4, TT14] . In our context, condition (2.14) 
is first needed to guarantee global-in-time well-posedness of solutions to (1.5), as discussed 
above. Then, condition (2.14) is needed again, in an essential way, to show that solutions fe 
of the kinetic model (1.3) remain compactly supported for all finite times, with a bound on 
their support that is independent of e (see Proposition 3.4 and its proof). This estimate is 
essential in the e ^ 0 limit of the kinetic model. 

Note that by [AGS06] (Lemma 8.1.6 in Chapter 8), the unique measure-valued solution 
in Theorem 2.4 is exactly the unique weak solution for (1.5), i.e., it satisfies 


J J x)p{t, x) dx dt 

0 Rd 

forany(^€Ci([0,r);Ci(]R‘')). 


T 

+ I I Vx 


ip-pudxdt-\- J (p{0,x)pq{x,v) dx = 0, ( 2 . 20 ) 


3 Well-posedness and uniform bounds of the kinetic model 

In this section we make preparations for proving the convergence of the kinetic equation (1.3) 
to the macroscopic model (1.5) by showing several uniform bounds in e. These include an 
uniform bound on the support of the density function /e, along with uniform bounds on its 
moments. A key step is to show that a certain “second-order fluctuation” term vanishes with 
e, which is the main ingredient that allows the passage from the second-order kinetic equation 
(1.3) to the hrst order macroscopic model (1.5). 

3.1 Well-posedness setup for the kinetic model (1.3) 

Model (1.3) with e > 0 fixed is a particular case of the general class of models studied in 
[CCRll]. For completeness, we present briefly the well-posedness setup from [CCRll], which 
is similar in fact to the measure transportation framework used for the macroscopic equation 
in Section 2. 

The solution space for equation (1.3) is (^([O, T], Rc(I^'^xl^'^)), where represents 

the space of probability measures with compact support in x R'^. The approach is to 
consider the characteristic equations associated with (1.3) and define measure-valued solutions 
/e to (1.3) as the mass-transport (or push-forward) of the initial density /o along the flow 
defined by characteristics. 
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Specifically, the characteristic equations for (1.3) are given by: 


dx 

du _ 1 

dt e 


\x -y\)iv - v*)feit,y,v*)dydv* + VxK * p^{t,x) , (3.1) 


(a:,u)|^^Q = (xo,uo) € supp/o. 

In short-hand notation, one can write the characteristic system (3.1) as 

dx 


dt 

dv 

dt 


= V, 


= 'H[fe]it,X,v), 


(3.2) 


(^>'^)L =0 = (®o,?^o) e supp/o, 


where T-Llfe] denotes the entire right-hand-side of the x-equation in (3.1). 

Suppose that is a hxed function and the vector held ?t[/e](t,x,x) is locally Lipschitz 
in X and x. Then, standard ODE theory provides the local existence of the how map 71^)^ j 
associated to (3.2): 


(Xo, Xo) (X, X), (X, X) = {Xeit),Veit)), 

where (xe(t), Xe(t)) is the unique solution of (3.2) that starts at (xo,xo). 

Similar to Dehnition 2.1, we consider measure-valued solutions of (1.3) in the mass- 
transportation sense [CCRll]: 

Definition 3.1. Let T > 0 be fixed. A measure-valued function fe € C([0, T]; Rc(I^'^ x IK'^)) 
is said to be a solution of the kinetic equation (1.3) with initial condition fo G x IK'^) if 

Mt) = rliif,]#fo, (3.3) 

for each t G [0, T]. 

The push-forward in (3.3) is in the mass transportation sense: 

I ({x,v)Mfix,v)dxdv = J C{rillf^]{X,V))MX,V)dXdV, (3.4) 

M2d ]R2d 

for all C G ^^(M^ x M'^). 

The results in [CCRll, Section 4] apply to the kinetic equation (1.3) with e > 0 hxed. 
The well-posedness result that can be inferred directly from there is the following. 

Theorem 3.1 (Well-posedness of the kinetic model (1.3) [CCRll]). Assume the following 
properties on the influence function and the interaction potential K: 

(j) and VK are locally Lipschitz, and |(/>(x)|, |VR'(x)| < C(1 -|- |x|) for all x G 
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for some C > 0. Take an initial measure fo G x K'^). Then there exists a unique 

solution /e G C([0, oo), x ^'^)) of (1.3), in the sense of Definition 3.1, whose support 

grows at a eontrolled rate. Specifieally, there exists an increasing funetion Re{T) sueh that 
for all T > 0, 

supp/e(t) C for all t G [0,r]. 

We also note that by results in [AGS06, CDF+11], given by Theorem 3.1 is also a weak 
solution of (1.3), i.e., it satisfies 


1 

j j aMt,.,v)Ut,.,vVV.AvM 


0 R2d 


+ / / • vfe dx du dt 

0 R2d 


J V.ui’■ F[fe]fedxdvdt+ J ^f{0,x,v)foix,v)dxdv = 0, 


(3.5) 


0 R2d 


for any 'll; G C^{[0,T)-,C^{R'^ x M"*)). 

Remark 3.2. An essential assumption made throughout our work is that the interaction 
potential K is symmetric. For such potentials, the kinetie model (1.3) conserves linear mo¬ 
mentum through time evolution. To see this, one ean take in the weak formulation (3.5) the 
test function if (t,x,v) = ifi{t)vif 2 {v), where if i G C^{0,T) andvif 2 {v) G C^(M'^), with 1^2 = 1 
on the set {u G : {x,v) G suppf^}. Due to the symmetry of K, fj^ 2 d F[ff:]f„dx dv = 0, and 
we obtain that 

d f 

— / vfe{t,x,v)dxdv = 0, in the distributional sense. 


Throughout the paper we assume, with no loss of generality, that the initial density fo has 
zero momentum, i.e., 

J vfo{x,v)dxdv = 0. (3.6) 

R2d 

Henee, the solutions fe provided by Theorem 3.1 satisfy: 


J vfeit,X,v)dxdv = 0, 

R2'^ 


(3.7) 


for all finite times. 

3.2 Uniform bound of the support 

We prove here that starting from a compactly supported probability density function fo, the 
solution ff:{t, ■, •) to (1.3) remains compactly supported in (x, v) for any finite time. Moreover, 
the size of supp/^ is independent of e. 

We first prove a rough bound on the first moment^ of /g. 

^We adopt here a common terminology used for kinetic equations and refer to moments of /e the moments 
with respect with the velocity variable v. 
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Proposition 3.3. Let be the solution to (1.3), as provided by Theorem 3.1. Suppose >0 
and that the initial data /o has a finite first moment. Then for each finite time T > 0 and 
each fixed e > 0, /e(t) has a finite first moment for any t € [0,T]. The bound of the first 
moment only depends on T, e, and K. In particular, it is independent of ([>. 


Proof. By symmetry of the integral we have 

R2d ]R4d 

+ y * fi) fe dx dr; 

R2d 

= \ f (Pi\x-y\)iv-v*)- fe{t,y,v*)fi{t,x,v)dxdydvdv* 


+ y ^1 • * fi) fe dx dn. 


Notice that 


Hence, 


This gives 


; ^ \ I I 1 I *1 

(v — V ) • T-; — ■;—; = n + n — 

' \V\ X* / ' ' ' ' 1^1 


> 0 . 


_d 

dt 


y 1^ • F[f^]fi dx d-u > y 1^ • {VxK * fi) fi dx dx. 

R2d R2d 

y |x|/e(t,x,x)dxdx <y ly • */,)/,dxdx < yiVjKlI^o 


Therefore, for t € [0, T], we have 

y |x|/e(t, X, x) dx dx < y |x|/o(x, x) dx dx H- 


1 ^ 2^-^11 L °° ■ 


Hence the first moment of is bounded for each e and T, and the bound is independent of 

fi. □ 

The following proposition shows that ffit,-,-) is compactly supported and the size of its 
support is independent of e. 

Proposition 3.4 (Uniform bound on the support of /g). Suppose K is symmetric, VxK E 
VU^’°°(M'^) and fi E VU^’°°(]R^) is positive, non-increasing, and satisfies the slow decay condi¬ 
tion (2.14). Let /e be the solution to (1.3) with an initial data /o E x M*^) that satisfies 

(3.6). Then for each finite time T, ffit,-,-) has a compact support for any t E [0,T]. The 
size of the support only depends on T, fi, and K; in particular, it is independent of e. 
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Proof. Consider the characteristic equations (3.1) associated with (1.3). Using the bound¬ 
edness of the first moment of on the time interval [0, T], as shown in Proposition 3.3, we 
immediately derive from (3.1) that for each fixed e and T, has a compact support in (x, v). 

Our goal is to show that the size of this compact support is independent of e. 

To indicate the dependence on the initial location, we use the notation (xe(t; xo,vo),Ve{t-, xq, uq)) 
for the characteristic path that originates at {xo,vo). Consider now the entire family of char¬ 
acteristic paths initialized at points (xo,no) inside the support of /q. Since Xe{t', xo,vo) and 
Vf:(t] xo,vo) are continuous in (xo,no) for each t G [0,T], the maxima of the supports of x^ 
and Ve are attained, and one can define the sizes S^it) and V^it) of these supports: 


Denote 


S'e(t) = max |xe(t; xo, no)|, U(t) = max |n£(t; xq, no)|. 
(3;o,i>o)6supp/o (xo,'Uo)Gsupp/o 


(3.8) 


Then 


Denote 


D,{t) = S'e(O) -F J |ne(T)| dr. 

0 

Se{t) < DS)- 

$s(t,x)= / (j){\x - y\)p^{t,y)dy. 


Let VM{t) be the velocity with the largest magnitude at time t, that is, 

\vMit)\ = Veit), 

and let T/(t) = f)i2De{t)). Then the monotonicity of cj) gives 

0 < yit) < (/>(|x - y\) for any (x, •), iv, •) G supp/e(t, •, •)■ 
The time evolution of Ve is given by 


dVe 1 VM 


dt e \vm\ 


\Xe - y\)ivM - V*)feit,y,V*) dy dx* -h VxK * Pe 


(3.9) 

(3.10) 

(3.11) 

(3.12) 


< - 


^i\vM\ J 4>i\xe-y\)peit,y)dy - J fUxe - y\)v* feit,y,v*) dy dv* 


(3.13) 


+ - \\VxK\\^^ = RHS. 
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By the momentum conservation (3.7), we have 


RHS = I \vM\^eixe) “ ' J (</>(ke “ 2/|) “ r]{t))v* fe{t, y, V*) dy dv* j + ^ 




< - 


- i\vM\<^e{xe) -\vm\ J { 4 >{\xe - y\) - v{t))fe{t,y,v*) dy dv* 

V R 2 '^ 

vit) 


+ - \\V,K 


II 


I 1 L° 

Using the estimate above in (3.13), along with (3.12) and the fact that rj{t) is non-increasing, 
we find 


d (y IIVx^IIlo 

dt V " vit) 


< (u, I 

< -I \Vm\ — 


< - 


e 

vit) 


vit) 


d / II Va;ii"|| 


dt \ r]{t) 


v,- 


l^a^-^ll L°° 


vit) 


Therefore for t € [0,7"] we have 


IV.KI 


Veit) < K(0) - " "".J, " + " "" - , < U,(0) + '' 


(3.14) 


7 ( 0 ) vit) vit) 

Combining (3.14) with the definition of D^it) in (3.9), we obtain the differential inequality 


Diit) < VM + ^4^1 = y^(0) + 


yit) ' <i>i2Deit)y 

Note that is an increasing function in t and (^(r) J- 0 as r —?> oo. To show that is 
bounded independent of e, we let 


r 

A(r) = J 


(pis) 


C2 + Cl(p{s 


■ ds, 


where ci = 14(0) and C 2 = ||Va;i7||^oo. Then A is an increasing function and it has the bounds 


C2 + Cl 


r 

- [ ^is)ds < A(r) 

I II<?>IIl°o j 


< — [ (j)is) ds. 
C2 J 
0 


By (2.14) we have Range(A) = [0,oo). Meanwhile, for each fixed e and T such that is 
well-defined on [0, T], solving (3.15) gives 

Ai2Deit)) < Ai2DeiO))+ 2t. 


Note that D^iO) = -S'e(O) is independent of e, as it represents the size of the spatial projection 
of the support of /o (see (3.8)). Therefore we obtain the following uniform in e bound: 

Deit) < ^A“^ (A(2|supp/o|) -b 2T) for any t G [0,r]. 

We conclude that by (3.10), we have that S^t) is uniformly bounded in e on [0, T]. Moreover, 
by (3.14) and the definition of y, we also have that the support of v is uniformly bounded in 
e on [0, Tj. □ 
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Throughout the rest of the paper we denote by ^{T) C x the common support of 
fe{t, •, •) for all e > 0 and t G [0,T]. Also, we use the following notation for the spatial and 
velocity projections of Q(T): 

Oi(r) = {x G : (x,u) G 0(r)}, f^2(T) = {u G : (x, u) G f^(r)} . (3.16) 


Using the uniform bound of v^, we can immediately infer the uniform (in e) bounds for 
all the moments of fe- 

Corollary 3.5. Suppose the assumptions in Proposition 3-4 hold. Then for any A; G N there 
exists a constant C{T,k,4>,K) such that 

/ \v\^ fe{t,x,v) dx dv < C{T,k,(l), K) < oo , k>l. (3-17) 

Remark 3.6. By doing a more careful estimate, one can show that (3.17) holds as long as the 
initial measure /o has bounded moments (with no requirement on /o having compact support). 
Nevertheless, we content ourselves with Corollary 3.5 since the main setting here requires that 
the initial data is compact supported. 

3.3 Uniform bound of the fluctuation 

In this part we show a uniform in e bound on the second order fluctuation term defined below 
by (3.20). This bound (and in fact the rate of vanishing of the second order fluctuation with 
e) is essential to show the convergence of [JabOO, FS14]. The estimates here are a priori. 
However, one can apply a density argument (see for example. Section 4 in [FS14]) to show 
that they hold for measure-valued solutions /^(t, •, •) G Pc (M'^ X M'^) at each t G [0,T]. 

By previous estimates, we have a priori uniform-in-e bounds for all the moments and the 
size of the support of over an arbitrary finite time interval [0,r]. Hence, for T > 0 fixed, 
there exists a constant C'o(T) that depends on 4> and K (but independent of e) such that, for 
k = 1,2, 

J |u|^/c(t, X, u) dx du + |supp/c(t, •, •)! < Co(r), for all t G [0, T]. (3.18) 

K2d 

Also, there exists a constant C'o(T) independent of e such that 

(j){\x — y\) > C(}{T) > 0, for any x,y G supp/c, and t G [0,T]. (3.19) 

Theorem 3.7. Suppose the interaction potential K and the influence function (f satisfy the 
assumptions in Proposition 3.4- Let fe be the solution to (1.3) starting from an initial data 
fo £ 'PcO^'^ X 1^^) that satisfies (3.6). Let be the second order fluctuation defined by 

Iflt)= j \F[flffldxdv. (3.20) 

R2d 

Then for any t G [0,T], there exists C{T) independent of e such that 

/c(t) <C(r)e2 + /c(0)e-^. 
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Proof. We first make the simple observation that by Cauchy-Schwartz and (3.18), for A: = 0,1, 
we have 


j \v\^\F[fe]\fedxdv< 
]R2d 



\v\‘^^f,dxdv 


1/2 

< Co{T) 


(3.21) 


Now we estimate the rate of change of 1^. To this end, we multiply ^ |T[/e]|^ to equation (1.3a) 
and integrate in x,v. This gives 

= \ j vfe-^x{\F[fe]\‘^'^ dxdv + ^ J (dt\F[ffi\'^'^ f^dxdv (3.22) 


Ri 


R2 


1 

Ye 


V„ 


]n-{F[m) dxdu. 


--V-- 

Do 

Below we estimate separately the terms Ri,R 2 and Dq indicated above. 

Estimate of Dq: We have 

Do = -^ J V„F : (F (8) F)f, dx du 

K2d 

= J - y)\F[fe]\‘^ fedxdv, (3.23) 

K2d 


where for the second equality we used (1.3b) and the definition (3.11) of d>e(x). By (3.19) we 
conclude that Dq provides dissipation in (3.22); the exact expression of Dq is used below to 
control a certain contribution from R 2 . 

Estimate of Ri: Erom the expression of Ri and (1.3b), we infer immediately: 

Ri< J |V,F[/,]||F[/,]||u|/,dxdi;, 

R2d 


where 

\'^xF[fe]\ < ||Va;(/>||j;^oo(-Rd) ^|u| + ||u/£||j;^l(-]g 2 d)^ + 11 11 . 

Hence, by (3.18) and Cauchy-Schwartz (see also (3.21)), there exists C'i(r) independent of e 
such that 

Ri<Ci{T)lF^{t). (3.24) 

Estimate of F 2 : The estimate of R 2 is more delicate. The expression of R 2 and (1.3b) yields 

R2= j {F[f,]-dtF[f,])f,dxdv 

R2d 

= F2I + R22, 
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where 


R 21 = j (l){\x -y\){v -v*)dtfe{t,y,v*) ■ F[f^]{t,x,v)fe{t,x,v)dxdydvdv* 

R4d 

and 

R22= j * dtPe{t,x) ■ F[f^]{t,x,v)feit,x,v)dxdv. 

We first estimate i? 22 - Integrating (1.3a) in v gives 

dtPe = -Vj; • (vfe). 

Angle brackets denote integration with respect to v. 

Hence by (3.18) and (3.21), there exists C' 2 (T) independent of e such that 

1^22! < j \V^K * {vfe)\ \F[f^]\{t,x,v)feit,x,v)dxdv 

< (3.25) 

Next we estimate i? 2 i- Using the kinetic equation (1.3a), we have 
R 21 = - J Hlx - y\){v - V*) {v* ■ Vyfe{t,y,v*)) • F[f^]{t,x,v)fe{t,x,v)dxdydvdv* 

R4d 

+ ^ y (t>{\x - y\){v - v*)Vy* ■ (U[/e]/e) {t,y,v*) • F[f^]{t,x,v)fe{t,x,v)dxdydvdv* 
R4rf 

= -^ 21,1 + -^ 21,2 • 

Integration by parts in i? 2 i,i gives 

i? 2 i,i = J Vy(/>(|x - y\) • {v*fe{t,y,v*)) (v - v*) ■ Flf^j(t, x, x, v) dx dy dv dv*. 

R4d 

Therefore, by (3.18) and (3.21), there exists CslT) independent of e such that 

ii?2i,ii < (ii^/.iIli(r2.) \\vF[fmLHR^^) + III II^[/.]/.IIl1(r2<^)) (3-26) 

<C3{T)lR^{t). 

We are left to estimate i? 2 i, 2 - Integration by parts in v* gives 

^ 21,2 = ~J 4>{\x - y\) {F[f^]f^) {t,y,v*) • (T[/,]/e) (t,x,n)dxdydndn*. 

R4d 
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Combining -^ 21,2 with Dq we get 

Do + R21,2 = \ j “ y\) ■ Ff*f^(lxdydv<lv*, 

R4d 

where F* = F[ff\{t,y,v*) and f* = ff:{t,y,v*). By symmetry, (3.19) and the conservation 
property /g 2 d .P"[/e]/e dxdr; = 0, we can infer further: 

-Do + ^^ 21,2 = -^ J (l){\x - y\)\F* - Fl"^ f*f^dxdydvdv* (3.27) 

K4d 

<-lCo(r) j \F*-F\^f:f,dxdydvdv* 

K4d 

= -^Co(r) j \F[f,ff,dxdv. 

K2d 

Putting together (3.24), (3.25), (3.26), and (3.27) we get 

Ri + R 2 + Do< C 5 {T)lF\t) - ^Co{T)R{t). (3.28) 

The following differential inequality for R can now be obtained from (3.22) and (3.28): 

Ipe < C,{T)iy\t) - ^Co{T)R{t), for t G [0,r]. 

By Gronwall’s inequality, we thus have 

R{t) < C{T) + 7,(0) e--^\ 

for some constant C'(T) independent of e. □ 

4 Passage to the limit in the kinetic equation 

In this section we show the convergence of solutions of the kinetic equation (1.3) to solutions 
Pe of the macroscopic equation (1.5), as well as the convergence of the characteristics of (1.3) 
to those of (1.5). 

4.1 Convergence to the macroscopic equation 

The convergence to the macroscopic solution is given in the following theorem. 

Theorem 4.1 (Convergence to the macroscopic model). Suppose K is symmetric, \/xK G 
1T^’°°(M'^) and (j) G 1T^’°°(]R''') is positive, non-increasing, and satisfies the slow decay con¬ 
dition (2.14). For any fixed T > 0 and e > 0, let G C([0, T];Pc(l^'^ x 1^“^)) be the unique 
solution to (1.3) with initial data fo G VcO^'^ x K"^) that satisfies (3.6). Then there exists a 
unique p G (^([O,T];Pc(M'^)) such that 

pfit, •) .) in p(]R'^) as g 0 , for each t G [0, T). 

Moreover, the limiting measure p is the unique solution to (1.5) with initial density po = 

/Rd/odr. 
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Proof. We apply a similar argument as in [JabOO] (see [FS14] as well). The key ingredient is 
to study the limit involving the term ^[/g] defined in (1.3b). For ease of notation, denote 

Je{t,x) = J vfeit,X,v)dv . (4.1) 

R'^ 


Then F[f^] can be written as 

F[f^]{t,X,v) =V(f*Pe{t,x) - (f* Je{t,x) + S/xK * Pe{t,x). (4.2) 

We divide the proof into four steps. 

Step 1. First we study the passage to the limit of J^. To this end, take the test function 
in (3.5) as il:{t,x,v) = 'ipi{t)'tp 2 {x) V'if 3 {v) where V’l G C^{0,T), 'tf 2 {x),v'if 3 {v) E Cj(M'^), and 
V'a = 1 on n 2 (T) (recall that ^ 2 iT) is defined in (3.16)). Using this if in (3.5), we have 


T 

j '*/’i(0 / lp 2 {x) Je{t,x) dx dt + 
0 Rd 


T 


[ Mt) [ 

0 R2d 


'^x'4’2 ■ V (8* vfe dx du dt 


T 

f ifift) f if 2 {x)F[f^]f^dxdvdt = 0. 
0 R2d 


Therefore, if we denote 



lf 2 {x) Je{t,x) dx. 


then its weak derivative is given by 

dt) = J '^x'if2-v®vfedxdv-^ j if2{x)F[f^]f^dxdv. 

R2d R2d 

By Theorem 3.7, there exists a constant C{ti,T) independent of e such that 

i I \F[f,]\f,dxdv<-^iUt)f^<Cih,T), 

R2d 

for any t E [ti, T] where 0 < ti < T. Using this fact, together with Corollary 3.5, we conclude 
that E kF^’°°(ti,T) for any V ’2 € and 

||6llu/i'<-(ti,r) ^ C{ti,T) ||'02|lrvi,oo for any ti > 0. 

Hence, {?e(i)}£>o is uniformly bounded and equicontinuous in C'([ti,T]), and by Ascoli-Arzela 
Theorem it converges uniformly on a subsequence. We conclude that for any if 2 E C'^(M'^) 
and ti E (0,T], there exists a subsequence e*, and ^(t) such that 

j if 2 {x)Jek{t^x)dx ^ i{t) inC'([ti,r]) as 0. (4.3) 

Rd 
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On the other hand, Proposition 3.4 provides a uniform (in e) bound for the support 
of fe, which in turn gives a uniform bound of Hence, for each t E [0,T), there exists 
J{t, •) E A4(M'^) and a subsequence of denoted as (t, •), such that 

^ in as Cfc;0. 

Note that in general, e^i could depend on t. However, comparing (4.3) with (4.4), 
that along the subsequence which is independent of t, 

J x) dx ^ J ip 2 ix)J{t,x) dx inC([ti,r]) as —>■ 0. (4.5) 

Rd, Rd 

By a density argument applied to we then have 

^ in7W(M'=') asefc^O, (4.6) 

where t E [ti,T], ti > 0 is arbitrary, and is independent of t. 

Furthermore, the assumption cf) E VF^’°°(]R‘^) guarantees that the sequence {(j)* Je}e>o 
uniformly bounded and equicontinuous in C{[ti,T] x Hi(T)), where Hi(T) is defined in (3.16). 
We conclude that along a subsequence which is still denoted as e^, 

(f) * Jef. ^ (f) * J in C ^[ti, r] X Hi(r)^ as Cfc ^ 0, (4.7) 

where again, > 0 is arbitrary, and is independent of t. 

Step 2. One can use similar arguments to infer convergence results for the other two terms 
entering the expression of F[f^] in (4.2). These arguments would in fact be identical to those 
used in [FS14] for the attractive-repulsive model with no alignment. We summarize the key 
ideas briefly and refer to the proof of Theorem 5.1 in [FS14] for details. 

By taking a test function of the form 'ip{t,x,v) = 'tpi{t)'ij) 2 {x) in (3.5), one can show, by 
similar calculations as in Step 1, that J^aip 2 {x)pe{t,x) is uniformly bounded in 
Together with the tightness of the sequence Pe{t, •) and the application of Prokhorov’s theorem 
(cf. [Bil71, Theorem 4.1]), it can be inferred that there exists a probability measure p{t, •) E 
P(M'’*) such that, along a subsequence of (still denoted as e^), 

Pe^{t,-) ^ p{t,-) inP(R‘^) as efc ^ 0. (4.8) 

The convergence above holds for each t E (0, T] and the sequence is independent of t. 

Also as in Step 1, one can get 

(j) * Pe^ ^ (j) * p in C X Hi(r)^ as Cfc ^ 0, (4.9) 

where the uniformity with respect to x follows from the equicontinuity in x of {0 * /?e(t, •)}. 

Finally, regarding the convergence of {^xK * Pej.}, the equicontinuity in x is immediate 
(by the assumption made on the potential). The equicontinuity with respect to time is 
slightly more delicate, as does not have enough regularity to be used as a test function. 
Nevertheless, one can overcome this difficulty by a regularization of the kernel (as shown in 
[FS14]) and it can be inferred that 

VxK * p^^^-^VxK * p in C X Hi(r)^ as ^ 0, (4.10) 


(4.4) 
we have 
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with ti > 0 arbitrary. 

Step 3. We now proceed to deriving the limiting equation (1.5b) at each t G [0,T). For all 
{x,v) G ll(r), we use (4.2), (4.7), (4.9) and (4.10) to derive that 

•) •, •) strongly in Ch{^{T)) for each t G (0,r], 


where 


Fo(t, x,v) = V (j) * p{t, x) — (j) * J{t,x) + VxK * p{t, x). 
Therefore, for each time t G (0,T], 

/ im,] - Fof fe, dxdu < ||F[/,J - ^0 as Cfc ^ 0. 

R2d 


Combining the above convergence with the vanishing in e of /^(t), as established in Theo¬ 
rem 3.7, we have 

j Fq dx du ^ 0 as 0. 

R2d 

Let Pi G C'^(M'^) be arbitrary. Then, by Cauchy-Schwarz inequality, 

J pi{x)FQf^^dvdx ^ 0 asek^O, 


and hence. 


J Pi{x) (4>Je^ - Pe^{(p* J + VxK * p)) dx 0 as ek^ 0. 


where ^ = (j) * p £ C'^(M‘^). By (4.6) and (4.8) we can then pass to the limit and get 
J pi{x) (<hJ — p{4>* J + VxK * p)) dx = 0 for any pi G (M*^). 


This implies that for each t G (0,T], 


4>(t,x)J(t,x) = p(t,x) ( j (t)(\x-y\)J{t,y)dy-VxK * p{t,x)\ . (4.11) 

i 


By denoting 

u{t,x) = J{t,x)/p{t,x), on the support of p(t, •), (4-12) 

then equation (1.5b) holds on suppp. 

Step 4- In this last step we show that the limit p satisfies the macroscopic equation (1.5a), 
as well as the momentum conservation (1.5c). First, following arguments used in [FS14], one 
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can generalize (4.5) and show that ■02 can be allowed to depend on t as well. Specifically, the 
following convergence holds: 

J a:) dx ^ J 'ip 2 {t,x)J{t,x) dx inC([ti,r]) as —>■ 0, (4.13) 

R'* R'^ 

for all e Cc{[ti,T),Cl{R'^)). 

Similarly, an analogous result can be shown for the convergence of 

j'il) 2 {t,x)pe^{t,x)dx ^ j 'il’ 2 it,x)p{t,x)dx in C {[ti,T]) as €k ^ 0, (4.14) 


for all V ’2 GC,([fi,r),Ci(K'^)). 

Now, in the weak formulation (3.5) for /^j., take a function ^ of the form ^(t,x,u) = 
(p{t,x) G Cc^([0,r),C'j((R‘^)) to get 

rr. J, 

/ (p(0, x)po(3:, v) dx = 0, (4.15) 

0 Rd 


0 Rd 


where po = fo(x, v) du. 

Fix ei > 0 small. Then, one can choose ti > 0 small enough such that 


ti ti 

J J \dt(pit,x)pe^{t,x)\ dxdt + J j dxdt<C 


w, 


1,00 ti < ei. 


(4.16) 


With ti fixed, such that (4.16) is satisfied, break the time integrals in (4.15) into 
two pieces: and f^. Due to (4.13) and (4.14) one can pass to the limit ^ 0 in the 

time integrals Hence, we infer from (4.15) and (4.16) that for every ei > 0, there exists 

0 < ti < ei/C\\p\\^i ,oo such that 


J J dt(p{t,x)p{t,x)dx dt + J J VxP ■ J dx dt + J {p{0,x)po{x,v) dx 


*1 R"* 
We conclude 
T 


< ei- 


11 dMt,x)p{t,x)dxdt + (p * J dx dt 4“ 9^(6, x^pq(^x^ u) dx — 0, (4.17) 


0 Rrf 


0 Rd 


for all (p{t,x) G ([0, T), C'f((M'^)). Note that (4.17) is equivalent to (2.20) since J = pu on 
the support of p and J = 0 outside the support of p. 

The conservation of momentum (1.5c) holds since it holds for every Indeed, by (3.7), 
we then have 

/ Je^(t,x)dx = 0 for all t G [0,T] and > 0. 

Rd 
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Consequently, by (4.5) and (4.12), we infer that (1.5c) is satisfied. 

We conclude by noting that the convergence above was derived on a subsequence 
However, by the uniqueness of the weak solution to (1.5), as established in Theorem 2.4, the 
full sequence pe(t, •) converges to p{t, •) for each t G [0,T]. □ 


4.2 Convergence of characteristic paths 

In this part we investigate in more detail the convergence of f^. In particular, we show that 
the characteristic paths of (1.3) converge to those of (1.5). This provides a geometric point 
of view of the singular limiting process, which also provides an explicit formula for the limit 
of fe- 

The characteristic system (3.1) along which the solution is transported has the form 


dx 

du 

e— = -(/)* p^{t, x) v + f)* Je{t, x) - VxK * pe{t, x), (4.18) 

(^>^)|i=o = (®o,^'o) G supp/o. 


The goal here is to pass the limit e —)■ 0 in the characteristic system (4.18) and relate the 
limit to the characteristic paths of the limiting macroscopic equation (1.5). 

The main tool in showing this limit is a classical result in singular perturbation theory 
due to Tikhonov [Tik52]. A brief account of this result is the following. Using notations from 
[Vas63], consider the general system 






T{x,v,t), 


(4.19) 


where x, u G and e > 0. 

System (4.19) is a two-scale equation, where t and r = t/e represent the slow and fast 
time scales, respectively. The slow dynamics of (4.19) is given by 


' dx 

, V = T{x,t), 

where v = T{x,t) is a root of the equation 

iF{x, V, t) = 0. 


(4.20) 


(4.21) 


Roots V = T{x,t) of (4.21) are in general non-unique. For a fixed configuration x* and time 
t*, the fast dynamics is defined by 

(Ill 

-=F{x\v,e). (4.22) 

Note that in the fast system (4.22), x* and t* are being regarded as parameters. 
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Tikhonov’s theorem establishes conditions on a root v = T(x,t) which guarantee that 
solutions of the two-scale system (4.19) converge, via a fast initial layer governed by (4.22), 
to solutions of the degenerate/slow system (4.20). 

A root V = T{x,t), with T defined on a closed and bounded set D C is called isolated 

if there is a 5 > 0 such that for all {x,t) G D, the only element in B(T{x,t),6) that satisfies 
B{x, v,t) = 0 is V = r(x, t). An isolated root T is called positively stable in D, if v* = T{x*,t*) 
is an asymptotically stable stationary point of (4.22) as r —)• oo, for each {x*,t*) G D. The 
domain of influence of an isolated positively stable root T is the set of points {x*,v,t*) such 
that the solution of (4.22) satisfying u|t-=o = v tends to v* = T{x*,t*) as r — oo. 

Tikhonov’s theorem [Tik52] states the following: 

Theorem 4.2 (Tikhonov [Tik52, Vas63]). Assume that a root v = T{x,t) of (4.21) is isolated 
positively stable in some bounded closed domain D. Consider a point (xq, uq, to) the domain 
of influence of this root, and assume that the slow equation (4.20) has a solution x{t) initialized 
at x{to) = xq, such that {x{t),t) lies in D for all t G [to,T]. Then, as e —)• 0, the solution 
{Xf:{t),Vf:{t)) of (4.19) initialized at {xo,vo), converges to {x{t),v{t)) := {x{t),T{x{t),t)) in 
the following sense: 

i) limue(t) = vit) for all t G (to,T*], and 
£->■0 

a) limxe(t) = x{t) for all t G [to,T*], 

£—>■0 

for some T* < T. 

Remark 4.3. The convergence of Ve{t) to v{t) occurs via a fast initial layer and normally 
does not occur at the initial time to, unless the initial data satisfies vq = r(xo,to)- 

Our main result concerning convergence of characteristic paths is the following theorem. 

Theorem 4.4 (Convergence of characteristic paths). Assume K and satisfy the same 
assumptions as in Theorem 4-1. Consider the measure-valued solution f^ to (1.3) and the 
characteristic path {xfit), vfit)) that originates from some (xo,uo) G supp/o at t = 0. Then, 


lim(xe(t), Ue(t)) = (x(t), u(t, x(t))) for all 0<t<T, (4.23) 


where x{t) is the characteristic trajectory of the limiting macroscopic equation (1.5) that starts 
at Xq, andu{t,x) is the velocity field defined in (4.12). In particular, x{t) satisfies (2.1), and 
u solves (1.5b). 

Proof. The characteristic system (4.18) has a right-hand-side that depends on e and Tikhonov’s 
theorem does not apply directly. To circumvent this, we replace {p^, Jfi by {p, J) in (4.18) to 
arrive at the following system: 


dx 

du 

e-^ = —(f) * p{t, x) V -\- (f>* J{t,x) — VxK * p{t, x), (4.24) 

(^’^)li=o = i^o,vo) G supp/o. 


Theorem 4.2 applies to system (4.24). Indeed, (4.24) can be written in the form (4.19), 
with 


T'(x, V, t) = —(f) * p{t, x) V + (f>* J{t, x) — VxK * p{t, x). 
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For a fixed spatial configuration x and time t, the root v = r(x,f) given by 


= A —71 —- V^K * p{t,x)) (4.25) 

is unique, hence isolated. It is also immediate that for a fixed spatial configuration x* and time 
t*, the corresponding fast equation (4.22) has a globally attracting equilibrium v* = T{x*,t*). 
Consequently, v* is positively stable and its domain of influence is {x*} x x {t*}. 

Denote by {xe{t),Ve(t)) the solution of (4.24) that originates from (xo,xo). Then, by 
Theorem 4.2, the convergence in (4.23), which needs to be shown for Xe{t) and Ve{t), holds 
for Xeit) and Xe(f) (note that by (4.11), the root T in (4.25) is in fact the velocity u defined 
by (4.12)). Hence, it would be enough to show that for a fixed f > 0, 

lim \xAt) — Xe(t)\ = 0 and lim \ve(t) — Ue(t)| = 0. (4.26) 

Indeed, from (4.18) and (4.24) we get 

e-^{Ve{t) - Ve{t)) = -(p* p^{t, Xeit))Veit) + (p * p{t, Xe{t))Ve{t) 

+ (p* J^{t,Xe{t)) -p* J{t,Xe{t)) - VxK * p^{t,Xe{t)) + VxK * p{t,Xe{t)). 

To the first line on the right-hand-side we add and subtract p* pe{t, x^{t))ve{t). Then, we can 
write 

e-^(^^e(i) - %{t)) = -p* Pe{t, Xeit)) (n,(f) - Veit)) + Geit), (4.27) 

where 

Geit) = ip* Peit, Xeit)) -p* pit, Xeit))) Veit) (4.28) 

+ p* Jeit, Xeit)) - p* Jit, Xeit)) - V xK * Peit, Xeit)) + V xK =f= p(t, Xeit)). 

By integrating (4.27) one finds 

t 

Veit) - Veit) = P*MTMr))dr ^11 t G [0, T]. (4.29) 

0 

Using Proposition 3.4, we have 

p * PeiT,XeiT)) > Cl, for all T G [0, T], 

where the constant Ci depends on T, p and K, but not on e. Consequently, from (4.29), we 
get 

t 

\veit) — Veit)\ — ~ J ^ |G'e(s)|ds, for all t G [0,7]. (4.30) 

0 

We now focus on estimating the right-hand-side of (4.30). Inspect Ge given by (4.28), 
in particular the term on the first line of the right-hand-side. The term Veit) is uniformly 
bounded in e (to show this, one can follow for instance the arguments used to prove Proposition 
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3.4). To estimate the term in round brackets, add and subtract Xe(t)). Then, by triangle 

inequality and the Mean Value Theorem, we get 

\(l)* Pe{t,Xe{t)) - (f)* p{t,Xe{t))\ < sup \\4> * {Pe - p)\\{ni{T)) 

te[o,T] 

+ sup ||V0*/9||£,oo(Q^(r))|Xe(t) - Xe(t)| 
te[o,T] 

Estimates entirely similar to the one above can be made for the terms on the second line of 
the right-hand-side of (4.28), which lead to: 

\4>* Je{t,Xeit)) - (p* J{t,X^{t))\ < sup \\(l) * (Je - J)\\L^{nr{T)) 

tG[0,T] 

-h sup \\V(t)* J\\L^(^n^(^T))\Xe{t) - Xe{t)\, 
tg[0,T] 

and 


VxK * p^(t,Xeit)) -V^K * p{t,Xe{t))\ < sup || * (Pe “ p) ||Lo°(ni(T)) 

4G[0,T] 

+ sup \\VlK * p\\Loo^Q^^T))\x^{t) - X^{t)\. 
t&[0,T] 


where fli(r) is dehned in (3.16). 

By the uniform convergences of * pe, (j) * J^, and VxK * Pe, as established in the proof 
of Theorem 4.1, and by assumptions we made on (j) and we can group the three estimates 
from above and get from (4.28), 


|G,(t)| <C, + C 2 |x,(t)-x,(t)|, 


(4.31) 


where Cg and C 2 are constants that depends on T, (p and K. Moreover, ^ 0 as e —>■ 0. 
Applying (4.31) in (4.30), we get 


t 

\ve{t) — v^{t)\ < ^ + — fe ^*'*“^^|xe(s) — Xe(s)|ds, for alH € [0,T]. 
El e J 


Using (4.18) and (4.24), we further find 

t S 

|xe(t) — Ve{t)\ < ^ + — J e J |ue(r) — fo(T)| dr ds, for all t £ [0,T]. (4.32) 

0 0 

Change order of integration in the double integral on the right-hand-side of (4.32) to get 

t t 

\ve{t) — Ve{t)\ <^+ — j\ve{T) — Ve{T)\ Js ds dx, for all t G [0, T]. 

0 r 

After evaluating the integral in s we finally arrive at 


t 

Ve{t) - fo(t)| <^ + ^J\v,{t) - fo(T)|dr, 

0 


for all t £ [0, T]. 


The second limit in (4.26) now follows from the integral form of Gronwall’s inequality, given 
that Ce —0 as e —?• 0. The convergence of trajectories follows from here as well. □ 
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The convergence of characteristic paths yields the limiting flow map T* : xq —>■ x{t). It 
is convenient in the calculations below to use the notation x(t;xo) to denote the limiting 
characteristic path x{t) that starts at xq. 

The next result characterizes the limiting densities. 

Theorem 4.5 (Characterization of the limiting densities). The limiting macroscopic density 
p identified in Theorem 4-i is the push-forward of the initial density po ^2/ limiting flow 
map T*', 

P = r^#Po- (4.33) 

In addition, for each t G [0,T), converges weak-* to a probability density f{t, •, •) G x 

M'^); 

fe^ f in V{R^ X R^) ase^O. (4.34) 

The limiting density f, with first marginal p, is given explicitly by: 

f{t, X, v) = p{t, x)6{v - u(t, x)), (4.35) 

where {p,u) is the unique solution of (1.5). 

Proof. The first part, expressed by equation (4.33), follows from considerations made in The¬ 
orem 4.1. However, it also follows directly, as a consequence of the argument below. 

The limiting behaviour of was not explicitly stated or needed in Theorem 4.1, but 
follows by arguments similar to those used for and pe in the proof of Theorem 4.1. Let us 
sketch this argument briefly. 

Fix V’l G (7^(0,T) and 'tl >2 G C'^(M‘^ x M'^), and let (p{t,x,v) = ‘fii{t)'tjj 2 {x,v) in (3.5). Find 


/ V ’2 (x, x)/e (t, X, x) dx dx dt 


0 R2d 


[Mt) [ 

0 R2d 


Va;'02 • V - -V„V^2 • Fife] 


fedxdvdt. (4.36) 


Denoting by 


ie{t)= I 'if 2 {x,v)ffit,x,v)dxdv, 


then, by (4.36), the weak derivative of is given by 


L{t) = 


VxV'2 • V - -V„V’2 • F[f^] 


fidxdv G L°°(0,r). 


The right-hand-side of the equation above is bounded (the boundedness of the term that 
contains e follows from Theorem 3.7). 

Since is uniformly bounded in 1T^’°°(0,T), it converges uniformly on a subsequence. 
On the other hand, similar to the arguments used for and in the proof of Theorem 
4.1, we note that the sequence ffit,-,-) G V{R'^ x M'^) is tight, and hence, for each t G 
[0, T), ffit, •, •) converges weak-* as measures, on a subsequence e^, to a probability measure 
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f{t, •, •) G P(R'^ X M'^). Note also that the subsequence does not depend on t, due to the 
equicontinuity of ie{t) derived above. 

Hence, 

fe.it,;-) ^ in X R'^') as ^ 0, (4.37) 

which proves the convergence (4.34) on a subsequence. To show its convergence on the full 
sequence e ^ 0 we use the uniqueness of /, as derived from the arguments below. 

Since /^(t) = ]#/o by Definition 3.1, (3.4) holds for /g^,: 

J ax,v)f,,it,x,v)dxdv = J C{rllll/X,V))fo{X,V)dXdV, (4.38) 

R2d R2d 


for all C G C6(R'^ x R'^). 

By the weak-* convergence of the left-hand-side of (4.38) converges as —)> 0: 

J ({x,v)fei^{t, x,v)dxdv ^ J Cix,v)f{t,x,v)dxdv. 

R2d ]K2d 

Due to convergence of trajectories (4.23), the right-hand-side of (4.38) converges by Lebesgue’s 
dominated convergence theorem, 

J C{T^ll^^iX,V))MX,V)dXdV ^ J ax{t-X),u{t,xif,X)))MX,V)dXdV, 

R2d R2d 

as Cfc —>■ 0. Combining these two, we find 

J ({x,v)f{t,x,v)dxdv= J C{x{t-,X),u{t,x{t;X)))fQ{X,V)dXdV, (4.39) 

]jj2d ^2^ 

for all C G Cb{R^ X R'^). 

First note that (4.33) can be derived from (4.39). Indeed, choose C(x,v) = (p(x) in (4.39) 
to find 

J (p{x)p{t,x)dx = J ip{x{t; X))pQ{X)dX, 

Rd- R'^ 

for all if G C'fe(R‘^); this represents exactly the mass transport given by (4.33). 

Now, observe that (4.35) is equivalent to 

J f{t,x,v)C{x,v)dxdv = J C{x,u{t,x))p{t,x)dx, 

R2d Rd 

for all test functions C G C'fc(R'^ x R'^), which can be inferred immediately from (4.33) and 
(4.39). 

The unique explicit representation of the limiting density / implies that the convergence 
in (4.37) holds on the full sequence /e, as stated in (4.34). □ 
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5 Numerical implementation and large time behaviour 


In this section, we present and implement a numerical scheme for the macroscopic system 
(1.5). The numerical approach taken here follows closely the analytical considerations made 
in Section 2. In particular, it provides a discrete analogue for replacing the singular operator 
A by an invertible operator At. 


5.1 Discrete setting 

The numerical discretization of the evolution equation (1.5a) for the density p can be done by 
standard methods. For instance, spatial discretization can be performed by a finite volume 
method [CCH15] or by a semi-Lagrangian particle method [DHPRF04]. For this reason, we 
focus in this section on the numerical solution of the velocity equation (1.5b). The subtlety 
is that for a given p, this equation does not have a unique solution for u, and hence, a direct 
discretization would lead to a singular system. 

We illustrate this degeneracy in one dimension. To this end, let h he a fixed mesh size 
and Xi = ih be equally distributed nodes. The range of i is —N,... ,N, with N G N large 
enough such that the support of the density is within the computational domain [—xn,xn]- 
The time dependence is irrelevant for the discretization of (1.5b) and we drop it in the 
calculations below. 

Denote by pi and Ui, i = —N,... ,N, the numerical approximations of p{xi) and u{xi), 
respectively. Also, denote by p and u the column vectors containing these values: 

P = {P-N, ■ ■ ■, Pn)'^, u = (u-N, • • •, un)'^- 


We choose the midpoint rule to approximate the integrals in (1.5b); note however that the 
argument below can be adapted to apply to higher order quadrature rules. This gives 



y\)p{y)dy 


- k\h)pk, 

k 


j 4'i\xi-y\)piy)u{y)dy 

R 


- k\h)pkUk, 
k 


j K'{xi - y)p{y)dy 
R 


hY,K'{{i-k)h)pu. 

k 


Introduce notations 


cl)i = cj){\i\h), K[ = K'{ih), i = -N,...,N. 


The following symmetries, respectively antisymmetries, hold: 

cl)i = (l)-i, K' = -K'_i, for alli = -iV,...,A, 

where the antisymmetry of K'- follows from (1.4). Also, since the influence function cj) is 
assumed to be positive and non-increasing (in accord with hypotheses made throughout the 
paper), (j) has a lower bound on the computational domain and we have (jii ^ V: a.11 

i = where rj = <j){2x]\f)- Note that we made an abuse of notation here, as a 
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variable t], with a very similar meaning and role, was defined in (2.4); however we try to 
parallel the analytical considerations made in the proof of Proposition 2.1, and refrain from 
introducing unnecessary extra notations. 

The discretization of (1.5b) reduces then to solving the linear system 


Au = b, 


(5.1) 


where 


A = 


and 


Zlfc f>-N-kPk 



4>-n-{-n)P-n ■ 

■ (p-N-NPN 


Y.k 4’N-kPk_ 


_ 4’n-{-n)P-n 

fiN-NPN_ 


b= {b-N,---,bN)'^, hi =i = 


(5.2) 


(5.3) 


The matrix A is the discrete analogue of the operator A defined in (2.6). Similar to the 
continuous case, matrix A is singular, and solutions of the linear system are not unique. This 
fact is detailed in the following proposition. 

Proposition 5.1 (Existence and non-uniqueness). Matrix A given by (5.2) is singular with 
dim{Null{A)) = 1. Moreover, the discrete linear system (5.1)-(5.3) has infinitely many 
solutions. 

Proof. It is easy to check that the entries of each row of A add up to zero; therefore, A is 
singular. To show that dim{Null{A)) = 1 we simply note that if we remove the i-th row 
and the i-th column of A (for an index i such that pi 7 ^ 0), the remaining matrix is strictly 
diagonally dominant, and hence it has full rank. 

To prove existence of infinitely many solutions, as opposed to no solution, we need to 
check b E Range{A). Indeed, by symmetry of fii, p = 0, and hence 

NulfiPA) = span{p}. 

Finally, by the antisymmetry of K[, one has b A p, which, by Fredholm alternative, yields the 
conclusion. □ 

Using Propostion 5.1, we immediately have 

Proposition 5.2 (Uniqueness with momentum conservation). There exists a unique solution 
u of (5.1)-(5.3), subject to the (diserete) momentum conservation condition 


p'^u = 0 . 


(5.4) 


Proof. Although this result can be derived directly from Proposition 5.1, using the structure 
of Null{A) deduced in its proof, we show here a more constructive way, that is also more 
relevant to the numerical implementation. Let e = (1, • • • , 1)^. Then, a solution u of (5.1)- 

(5.3) which satisfies the discrete momentum conservation condition (5.4), also solves 


Mu = b, 


(5.5) 
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with 


(5.6) 


M = A + 7] ep^. 

Here, p represents the positive lower bound of (j)i, * = —Al,... ,N, mentioned above. 
The matrix M has the form 


Y^k 4'-N-kPk 


{4'-N-{-N) - V)p-N ■ ■ ■ ifi-N-N - P)PN 

Yk 4’N-kPk_ 


_ i^N-i-N) - P)P-N ■ ■ ■ {4>N-N - V)PN _ 


Since all the off-diagonal entries of M are negative, M can be shown to be strictly diagonally 
dominant. Hence, M is an invertible matrix, and (5.5) has a unique solution. □ 

Remark 5.3. The discrete theory aligns perfectly well with the results in Section 2. In 
particular, matrices A and M are the discrete analogues of the singular, respectively invertible, 
operators A and A4 defined in (2.6) and (2.7), while the conservation (5.4) of the discrete 
linear momentum is used as an auxiliary condition that enforces uniqueness of the numerical 
solution. 

To conclude, our numerical discretization of (1.5b) consists in the following procedure. 
First, form the matrix A using appropriate quadrature rules on integrations. Next, change 
the singular matrix A into the invertible matrix M given by (5.6), and solve the linear system 
(5.5). For an efficient numerical implementation, it is important to note that the velocities 
Ui at grid points with nontrivial density values pi 0 do not depend on values Uj for which 
Pj = 0; this comment is in fact the discrete counterpart of Remark 2.2. Hence, one can 
consider only the nodes (or indices) that correspond to nontrivial densities and reduce the 
size of the system. 


5.2 Numerical examples 

Below we illustrate our method with three numerical simulations: two examples in ID with 
smooth and non-smooth potentials, respectively, and a 2D example with a smooth potential. 
We make choices of potentials for which explicit steady states of (1.5) can be calculated. The 
influence function in all the cases is chosen as 

(?i(r) = (1-|-(5.7) 


Example 1 (ID with smooth potential). We consider a power-law interaction potential 
with quartic attraction and quadratic repulsion: 


K{x) 



(5.8) 


The derivative K' is smooth, but not bounded, which means that K’ lies in Wjjj’^(]R), but not 
in VF^’°°(M), as assumed by theoretical results in previous sections. In simulations however, 
we deal with a bounded computational domain, on which K' and K” are of course, bounded. 
The initial condition is chosen to be a smooth bump function 

Pq{x) = Cse~^^li^_s^s){x), (5.9) 
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where s is a positive integer which determines the support of poi 1 denotes the indicator 
function, and Cg is a normalization constant such that UpoIIli = 1 - 

Figure 5.1 shows the time evolution of p starting from initial data (5.9) with s = 0.5 
(top row) and s = 1.5 (bottom row). The computational domain is set to be [—2,2] and the 
mesh size is h = .02. In both cases we observe that the solution approaches (by expanding, 
respectively compressing) the same steady state poc given by two Dirac delta distributions 
located at —1/2 and 1/2: 

Poo{x) = ^ [(5(a: - 1/2) + 6{x + 1/2)] . (5.10) 

We note that (5.10) is in fact a steady state of the attraction-repulsion model (with no 
alignment), i.e, 

dtp-Vx ■ ip{VxK * p)) = 0. (5.11) 



Figure 5.1: Time evolution of the solution p{t,x) of (1.5) in one dimension, starting from the 
initial density (5.9) with s = 0.5 (top row) and s = 1.5 (bottom row). The influence function 
is given by (5.7) and the interaction potential by (5.8). The two density profiles expand, 
respectively compress, toward the same equilibrium state given by (5.10). 


Example 2 (ID with Morse potential). We now take a Morse-type interaction potential: 

K(x) =-e-l^l/2 + e-l^l. (5.12) 

Note again that due to the singularity at the origin, K' 0 1F^’°°(M). In terms of numerical 
simulations however, singularities in the potential are not an issue of concern (at discrete 
level, a pointy potential and a very localized regularization of it, are essentially the same). 

In this simulation we choose a more general non-smooth and non-symmetric initial density: 

/3o(a:) = ^(1 - x)l(_m)(x). (5.13) 
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Figure 5.2 shows the time evolution of the solution (solid line) at t = 0, 2,5, 20. In the 
same plot we show (dashed-line) the steady state of the attractive-repulsive equation (5.11) 
with the Morse-type potential (5.12). An exact, explicit form of this equilibrium solution was 
derived in [BTll]: 

Poo{x) = C {cos{p{x - c)) - A) • 1\x-c\<h{x), (5.14) 

where the constants are given as follows; p, = l/\/2, A = —1/(3 -|- H = Also, 

C = 3/(2(7r -|- 3\/2)) is a scaling that sets the mass to one, and c = 1/3 is a shift of the 
centre of mass to the origin. Similar to Example 1, we observe that the solution approaches 
asymptotically the steady state poo- 


Q, 


t = 0 


t = 2 


t = 5 


t = 20 






X X X X 


Figure 5.2: Time evolution of the solution p{t,x) of (1.5) in one dimension, starting from the 
initial density (5.13). The influence function is given by (5.7) and the interaction potential 
by (5.12). The solution (solid curve) approaches asymptotically the equilibrium state poo 
(dashed line) of the attractive-repulsive model (5.11), given explicitly by (5.14). 

Example 3 (2D with Newtonian-quadratic potential). In the final example, we check 
the asymptotic behaviour of solutions to (1.5) in 2D. The algorithm presented above can be 
immediately extended to two dimensions. We take an interaction potential K with Newtonian 
repulsion and quadratic attraction: 

iF(x) = — log |a:|-|-(5.15) 

With this choice of AT, it has been shown that the steady state of the attractive-repulsive 
model (5.11) is the indicator function of a disk: 

Poo(a;) = C'1 |,k|<i(x- c), (5.16) 

where the scale C is set by conservation of mass, and c is a shift of the centre of mass [FHKll]. 
We take the following radially symmetric initial density: 

Po(r) = (r - lf{r - 2fli^r<2{r). (5.17) 

The numerical time evolution illustrated in Figure 5.3 indicates again, that the solution 
approaches asymptotically poo (here C = 0.1 and c = 0). 

In this experiment we set the computational domain to be [—2.2, 2.2] x [—2.2, 2.2] and fix an 
Eulerian computational grid with dx = dp = 0.04. We use a semi-Lagrangian scheme to evolve 
the density equation in time [DHPRE04]. Specifically, we evolve numerically, using a forward 
Euler integrator with dt = 0.005, the characteristic trajectories (2.1) and the density equation 
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t=0 


t=0.6 


t=1.2 



Figure 5.3: Time evolution of the density p{t, x) of the aggregation model (1.5) in two dimen¬ 
sions, starting from the initial density (5.17). The influence function is given by (5.7) and 
the interaction potential by (5.15). The solution approaches asymptotically the equilibrium 
state (5.16) consisting of a constant distribution in the unit disk. 


(1.5a) written in characteristic form. After each time step, we interpolate the density values 
from the Lagrangian (distorted) grid back to the Eulerian grid, calculate the new velocity 
field according to the method described in Section 5.1, and repeat the procedure. 

We conclude this section by noting that our preliminary numerical investigations of model 
(1.5) (as presented here) suggest that the large time behaviour of solutions to (1.5) is: 

p{t,x) —)• pooix), for X E suppp{t), as t ^ oo, 

where poo is a steady state of (5.11). Consequently, 

u(t, x) —>■ 0, for all x E suppp(t), as t ^ oo. 

This means that the macroscopic velocity u converges to its (zero) average as time goes to 
inhnity, which corresponds to flocking. Whether the model always approaches asymptotically 
the flocking limit is left for future investigation. 


37 








References 


[AGS06] 

[AIRll] 

[BCL09] 

[BDF08] 

[Bil71] 

[BL07] 

[BLRll] 

[BTll] 

[BV06] 

[CCH14] 

[CCH15] 

[CCRll] 

[CD 10] 
[CDF+03] 

[CDF+11] 

[CDM+07] 

[CFRTIO] 

[CFTVIO] 

[CKJ+02] 

[CS07] 


Luigi Ambrosio, Nicola Gigli, and Giuseppe Savare, Gradient flows: in metric spaces and in the 
space of probability measures, Springer Science & Business Media, 2006. 

Martial Agueh, Reinhard Illner, and Ashlin Richardson, Analysis and simulations of a refined 
flocking and swarming model of Cucker-Smale type, Kinet. Relat. Models 4 (2011), no. 1, 1-16. 

Andrea L. Bertozzi, Jose A. Carrillo, and Thomas Laurent, Blow-up in multidimensional aggre¬ 
gation equations with mildly singular interaction kernels. Nonlinearity 22 (2009), no. 3, 683-710. 

Martin Burger and Marco Di Francesco, Large time behavior of nonlocal aggregation models with 
nonlinear diffusion, Netw. Heterog. Media 3 (2008), no. 4, 749-785. 

Patrick Billingsley, Weak convergence of measures: Applications in probability. Society for Indus¬ 
trial and Applied Mathematics, CBMS-NSF Regional Conference Series in Applied Mathematics, 
Philadelphia, Pennsylvania, 1971. 

Andrea L. Bertozzi and Thomas Laurent, Finite-time blow-up of solutions of an aggregation 
equation in R", Comm. Math. Phys. 274 (2007), no. 3, 717-735. 

Andrea L. Bertozzi, Thomas Laurent, and Jesus Rosado, theory for the multidimensional 
aggregation equation, Comm. Pur. Appl. Math. 64 (2011), no. 1, 45-83. 

Andrew J Bernoff and Chad M Topaz, A primer of swarm equilibria, SIAM Journal on Applied 
Dynamical Systems 10 (2011), no. 1, 212-250. 

M. Bodnar and J. J. L. Velazquez, An integro-differential equation arising as a limit of individual 
cell-based models, J. Differential Equations 222 (2006), no. 2, 341-380. 

J.A. Carrillo, Y.P. Choi, and M. Hauray, The derivation of swarming models: Mean-field limit 
and Wasserstein distances, Collective dynamics from bacteria to crowds, 2014, pp. 1-46. 

Jose A Carrillo, Alina Chertock, and Yanghong Huang, A finite-volume method for nonlinear 
nonlocal equations with a gradient flow structure. Communications in Computational Physics 17 
(2015), no. 01, 233-258. 

Jose A Canizo, Jose A Carrillo, and Jesiis Rosado, A well-posedness theory in measures for some 
kinetic models of collective motion. Mathematical Models and Methods in Applied Sciences 21 
(2011), no. 03, 515-539. 

Felipe Cucker and Jiu-Gang Dong, Avoiding collisions in flocks, IEEE Transactions on Automatic 
Control 55 (2010), no. 5, 1238-1243. 

Scott Camazine, Jean-Louis Deneubourg, Nigel R. Franks, James Sneyd, Guy Theraulaz, and Eric 
Bonabeau, Self-organization in biological systems, Princeton Studies in Complexity, Princeton 
University Press, Princeton, NJ, 2003. Reprint of the 2001 original. 

J. A. Carrillo, M. DiFrancesco, A. Figalli, T. Laurent, and D. Slepcev, Global-in-time weak measure 
solutions and finite-time aggregation for nonlocal interaction equations, Duke Math. J. 156 (2011), 
no. 2, 229-271. 

Yao-Li Chuang, Maria R. D’Orsogna, Daniel Marthaler, Andrea L. Bertozzi, and Lincoln S. 
Chayes, State transitions and the continuum limit for a 2D interacting, self-propelled particle 
system, Phys. D 232 (2007), no. 1, 33-47. 

Jose A Carrillo, Massimo Fornasier, Jesiis Rosado, and Giuseppe Toscani, Asymptotic flocking 
dynamics for the kinetic Gucker-Smale model, SIAM Journal on Mathematical Analysis 42 (2010), 
no. 1, 218-236. 

Jose A. Carrillo, Massimo Fornasier, Giuseppe Toscani, and Francesco Vecil, Particle, kinetic, 
and hydrodynamic models of swarming. Mathematical modeling of collective behavior in socio¬ 
economic and life sciences, 2010, pp. 297-336. 

1. D. Couzin, J. Krause, R. James, G.D. Ruxton, and N. R. Franks, Gollective memory and spatial 
sorting in animal groups, J. Theor. Biol. 218 (2002), 1-11. 

Felipe Cucker and Steve Smale, Emergent behavior in flocks, IEEE Transactions on Automatic 
Control 52 (2007), no. 5, 852-862. 


38 



[DCBC06] 

[DHPRF04] 

[FHKll] 

[FS14] 

[HL09] 

[HT08] 

[JabOO] 

[JaclO] 

[JE07] 

[KSUBll] 

[LLEK08] 

[LTB09] 

[MEK99] 

[MKJNR12] 

[MTll] 

[MT14] 

[Rey87] 

[SR14] 

[Tanl4] 

[TB04] 

[Tik52] 

[TT14] 

[Vas63] 


Maria R. D’Orsogna, Yao-Li Chuang, Andrea L. Bertozzi, and Lincoln S. Chayes, Self-propelled 
particles with soft-core interactions: patterns, stability and collapse, Phys. Rev. Lett. 96 (2006), 
no. 10, 104302. 

J. Donea, A. Huerta, J.-P. Ponthot, and A. Rodriguez-Ferran, Arbitrary Lagrangian-Eulerian 
Methos, Encyclopedia of Computational Mechanics, 2004, pp. 413-437. 

R. C. Fetecau, Y. Huang, and T. Kolokolnikov, Swarm dynamies and equilibria for a nonloeal 
aggregation model, Nonlinearity 24 (2011), no. 10, 2681-2716. 

Razvan Fetecau and Weiran Sun, First-order aggregation models and zero inertia limits, arXiv 
preprint arXiv:1410.7095 (2014). 

Seung-Yeal Ha and Jian-Guo Liu, A simple proof of the cucker-smale flocking dynamics and 
mean-field limit. Communications in Mathematical Sciences 7 (2009), no. 2, 297-325. 

Seung-Yeal Ha and Eitan Tadmor, From partiele to kinetic and hydrodynamic descriptions of 
flocking. Kinetic and Related Models 1 (2008), no. 3, 415-435. 

Pierre-Emmanuel Jabin, Macroscopic limit of Vlasov type equations with friction, Annales de 
rinstitut Henri Poincare (C) Non Linear Analysis 17 (2000), no. 5, 651-672. 

M. O. Jackson, Social and Economic Networks, Princeton University Press, 2010. 

M. Ji and M. Egerstedt, Distributed coordination control of multi-agent systems while preserving 
connectedness, IEEE Trans. Robot. 23 (2007), no. 4, 693-703. 

Theodore Kolokolnikov, Hui Sun, David Uminsky, and Andrea L. Bertozzi, A theory of complex 
patterns arising from 2D particle interactions, Phys. Rev. E, Rapid Communications 84 (2011), 
015203(R). 

Yue-Xian Li, Ryan Lukeman, and Leah Edelstein-Keshet, Minimal mechanisms for school forma¬ 
tion in self-propelled particles, Phys. D 237 (2008), no. 5, 699-720. 

Andrew J. Leverentz, Chad M. Topaz, and Andrew J. Bernoff, Asymptotic dynamics of attractive- 
repulsive swarms, SIAM J. Appl. Dyn. Syst. 8 (2009), no. 3, 880-908. 

A. Mogilner and L. Edelstein-Keshet, A non-local model for a swarm, J. Math. Biol. 38 (1999), 
534-570. 

Jennifer M. Miller, Allison Kolpas, Joao Plinio Juchem Neto, and Louis F. Rossi, A continuum 
three-zone model for swarms. Bull. Math. Biol. 74 (2012), no. 3, 536-561. 

Sebastien Motsch and Eitan Tadmor, A new model for self-organized dynamics and its flocking 
behavior. Journal of Statistical Physics 144 (2011), no. 5, 923-947. 

Sebastien Motsch and Eitan Tadmor, Heterophilious dynamics enhances eonsensus, SIAM Review 
56 (2014), no. 4, 577-621. 

Craig W Reynolds, Flocks, herds and schools: A distributed behavioral model, ACM Siggraph 
Computer Graphics 21 (1987), no. 4, 25-34. 

L. Saint-Raymond, A mathematical PDF perspective on the Chapman-Enskog expansion. Bull. 
Amer. Math. Soc. 51 (2014), 247-275. 

Changhui Tan, A discontinuous Galerkin method on kinetic flocking models, arXiv preprint 
arXiv:1409.5509 (2014). 

C. M. Topaz and A. L. Bertozzi, Swarming patterns in a two-dimensional kinematic model for 
biological groups, SIAM J. Appl. Math. 65 (2004), 152-174. 

A. N. Tikhonov, Systems of differential equations containing small parameters in the derivatives. 
Mat. Sb. (N.S.) 31(73) (1952), 575-586. 

Eitan Tadmor and Changhui Tan, Critical thresholds in flocking hydrodynamics with non-local 
alignment. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engi¬ 
neering Sciences 372 (2014), no. 2028, 20130401. 

A. B. Vasil'eva, Asymptotic behaviour of solutions of certain problems for ordinary non-linear 
differential equations with a small parameter multiplying the highest derivatives, Uspehi Mat. 
Nauk 18 (1963), no. 3 (111), 15-86. 


39 



